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Abstract 

A non-Hermitian A^— site-lattice Hamiltonian H with Laguerre-polynomial 
right eigenvectors and real energies is made self-adjoint in an ad hoc Hilbert 
space 1-L^^\ The necessary physical inner products are defined via alternative 
A;— parametric (2A; + 1)— diagonal metrics = 0^ constructed, in closed form, 
a:tk = Q,k = l,k = 2 and /c = 3. The value of k is interpreted as a degree of 
non-locality of the model. 



1 Introduction 



1.1 The Laguerre-polynomial quantum model 

It is well known that the Laguerre polynomialJ3 L(n, a, z) may formally be 
arranged in an infinite-dimensional Dirac-ket-like column vector 



|^{oo)^ 



L(0,a,^) (= 1) 
L(l, a,z) {= a + 1 — z) 
L(2, a, z) [= \{a + l){a + 2) - {a + 2) z + \ z^] 



(1) 



which satisfies, line-by-line and at any complex z G C, the infinite linear 
algebraic system of equations 



(2) 



which resembles the Schrodinger equation of quantum theory and where 



a+ 1 



-a-1 a + 3 -2 



-a — 2 a + 5 



-a — 3 a + 7 



(3) 



One can immediately formulate at least three reasons ag ainst i7(~)(a) of 
Eq. ([3]) being interpreted as a quantum Hamiltonian: 

• (pragmatic reason) the set of all complex z E C does not look like a 
good candidate for a phenomenologically meaningful spectrum; 



^our present denotation L{n,a,x) and normalization L{—l,a,x) — 0, L{0,a,x) — 1 
are taken from the symbolic-manipulation langauge MAPLE up to a factor the same 
polynomials are denoted as L'^{x) in Ref. |2] or as L'n\x) in Ref. [5] 
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• (mathematical reason) in the current, "friendly" Hilbert space £^ (to be 
denoted by symbol 'H^^\ cf. Appendix) the norm of the wave-function 
candidate ([T]) would be infinite; 

• (physical reason) the candidate H^°°\a) for the Hamiltonian of the 
system is manifestly non-Hermitian. 

1.2 The finite-dimensional Laguerre-polynomial quan- 
tum model 

Let us replace the matrix-resembling array H^°°^ (a) by its finite-dimensional, 
truncated version 



i/(^)(a) 



a + 1 -1 



-a — 1 a + 3 







-a — 2 a + 5 











-a-N + 1 a + 2N-l 



(4) 



A return of the finite-dimensional descendant of Eq. ([2]) 



//(^)(a) l^f )) = )(a) l^f )) , n = 1, 2, . . . , iV - 1 



(5) 



to the status of Schrodinger equation finds a new support in a few suddenly 
emerging encouragements: 

• (pragmatic encouragement) after truncation, the set of the admissible 
eigenvalues shrinks from the whole complex plane of z to the discrete 
and non-degenerate A^— plet of strictly real zeros of the elementary 
secular equation 



L{N,a,z) = 0, z = En, n = 0,1, . . . , N - 1 ; 



(6) 
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• (mathematical encouragement) after truncation (and in both the Hilbert 
spaces T-L^^'^^ defined in Appendix), the finite array 

L{0,a,En) 
L(l, a, En) 

L{N-l,a,En) _ 
may very easily be normalized, say, to one. 

An obstacle which seems to survive is the third, non-Hermiticity point and 
counterargument. In what follows we intend to show that this argument is 
misleading and that in a way summarized in Appendix the Hermiticity of 
Hamiltonian H^^\a) may be achieved and guaranteed in a large number of 
alternative physical Hilbert spaces H^^^ where the inner products are defined 
in terms of the so called metric operator 0. We shall also show that for 
our particular model these spaces (i.e., the metric-operator matrices) may be 
constructed by purely non-numerical means. 



Table 1: A sample of the dimension- and parameter-dependence of the energy 
spectra {i?„(a)} of Hamiltonian (jl]). 



parameters 


energies 




a 


Eoia) 


E,ia) 






En -lia) 


6 


1.0 


0.5276681217 


1.796299810 




11.23461043 


17.64596355 




2.0 


0.8899410156 


2.433144232 




12.60041387 


19.26204255 




3.0 


1.296419203 


3.093998381 




13.94134537 


20.83985455 


9 


1.0 


0.3681784529 


1.243357962 




20.38218199 


28.11834338 




2.0 


0.6318537723 


1.712163195 




21.90120660 


29.82533613 




3.0 


0.9343511232 


2.208578822 




23.39499254 


31.50012806 



l^fna)) = 



4 



1.3 A phenomenological addendum 



Our present solvable toy- model (jl]) is just a very special case of the broad 
class of the A^— dimensional tridiagonal-matrix Hamiltonians of the form 



Him 



ai 


Cl 











b2 


a2 













bs 


as 


cs 























bN~i 















bN ClN 



which describe the one- dimensional single-particle A^— site quantum-lattice 
dynamics reduced to the mere nearest-neighbor interaction. The kinematical 
aspect of virtually all of these models is, typically, characterized by the 
one-to-one correspondence between the matrix index s = 1,2, . . . , N and a 
discrete coordinate of the "site". In a more dimensional case these sites 
may be numbered by an ad hoc multiindex, s —> s, but just the simplest 
one-dimensional case is to be considered in what follows. 

At the very start of our present application of the theoretical innovation 
of the formalism using 7^ / in combination with Hamiltonian ([8]) and/or its 
special case (jl]) let us emphasize that such an effort may even be supported by 
a purely pragmatic motivation. The numerical sample of zeros z — En ^ (a) 
of Eq. as displayed in Table [T] persuades us, for example, that each of 
the A^— plets of the present toy-model bound-state energies really looks very 
much like a realistic phenomenological spectrum. 

Besides the similar intuitive argument it is equally important to keep in 
mind that the necessary modifications and adaptations or deformations of 
the spectrum are controlled not only by the freely variable coupling alias 
self-coupling alias asymmetry-strength a > but also by the very choice of 
the ad hoc model-space dimension A^ = 1, 2, . . .. 

This means that even on the level of practical applicability, our unusual 
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model (jl]) might find its place among the other standard fitting tools for the 
description of experimental quantum spectra, especially when they are not 
explained by the conventional interactions yielding vibrations (like harmonic 
oscillator) or rotations (mainly in more spatial dimensions). Thus, although 
the core of our present message should be predominantly methodical and 
mathematical, even its purely pragmatic aspects might prove relevant and/or 
useful. 



2 The concept of position in quantum lattice 



In the vast majority of phenomenological quantum-lattice models the Hamil- 
tonian is assumed Hermitian in the most common A^— dimensional vector 
space, i.e., we have H^^"' = (^H^'^^^ . In parallel, the coordinates Qs 
sumed real and observable so that, in the Schrodinger's "mode of description" 
[1] , the corresponding quantum operator of position is most often represented 
by a diagonal and time-independent A^— dimensional matrix 



\ 



Q2 



(9) 



Using the standard Dirac's notation the time-dependent Schrodinger equa- 
tion then reads 

\d, m)) = H^""^ \m) (10) 



and controls the evolution of the system prepared in a normalized pure state 
at time tprep = 0. At the time of measurement t > one evaluates the wave 
function {qsl'ipit)) = 'ip{t,s) and determines the probability 



g{t,s) = \{qMt))\^ = r{t,s)tP{t,s) 



(11) 



of finding the particle (or quasi-particle) at the s— th site g<j 
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2.1 The loss of the measurabiUty of position in the so 
called PT— symmetric Quantum Mechanics 

In the recent literature there emerged an interesting modification of the whole 
paradigm (cf. its compact summary in Appendix). It has been developed, 
mainly, within the framework of the so called PT— symmetric quantum me- 
chanics (cf., e.g., comprehensive reviews O [6] and also ^ or ^). In what 
follows we intend to describe one of applications of this new paradigm to 
our exactly solvable single-particle A^— site quantum-lattice model (jl]) where 
the spectrum remains real but where the Hamiltonian matrix itself appears 
manifestly non-Hermitian in i'^ = 'H^^\ 

The technical essence of the new paradigm may be found summarized in 
Appendix. In nuce, the key idea of the new formalism lies in the weakening 
of the standard textbook Hermiticity H = //^ of the Hamiltonian and in its 
replacement by the requirement 

H^Q = QH. (12) 

The use of this relation containing a nontrivial metric operator Q ^ I has 
already been shown useful, almost twenty years ago, in nuclear physics [7]. 
Naturally, the same requirement applies also to the other operators O of ob- 
servables where the rule = 0'^ must again be replaced by its generalization 

e o = ot e . 

In opposite direction, the "lattice-position" matrix ([9]) loses its status of 
an observable quantity whenever the metric ceases to be a diagonal matrix. 
This argument (nicely explained also in Ref . [H] ) ) further implies that also the 
quantity g{t, s) of Eq. (ITT]) will lose its original probabilistic interpretation 
in general. 
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2.2 A return to the (smeared) positions. 

A partial resolution of the puzzle of the smearing of positions was offered 
in Ref. [10] where the concept of the position and locality has partially been 
restored via the use of the "next to diagonal" , band-matrix, ( 2 A; +1)— diagonal 
matrices of the metrics B = Gfc = 

Oil 9i2 ... ^i,fc+i ... 

^21 ^22 ^23 • • • ^2,A:+2 ' • '■ 

\ ■•. ■•. ■•. ■-. 

dk+l,l ' ■ ' ■ ' ■ ON-k,N 

■ . ■ . On-2,N~2 0n^2,N~1 '■ 

■ ' ■ (^N-l,N-k-l ■ ■ ■ 0N-l,N-2 ^Af-l,Af-l ^Af-l,Af 

... ON,N~k ■ ■ ■ 9n^n~i 6nn 

"(13) 

The integer A; = 0,1,...,A^ — 1 has been interpreted there as a measure of 
the "smearing" . In this role, its value must be kept perceivably smaller than 
the dimension A^. 

The physical meaning of the above non-diagonal metrics may be clarified 
using the first nontrivial k = 1 example in which we may rewrite metric (fT3|) 
as a superposition of a positive diagonal matrix with and upper-diagonal 
matrix a A and its lower- diagonal conjugate aA\ 

e[^^ (a, a) = V^ + a {A + A^) . (14) 

Once we assume that the parameter a itself is sufficiently small, we may recall 
formula fHT]) of our Appendix and deduce, say, the following approximate, 
non-Hermitian form of the Dyson's map. 



n = ^e[''\a,a) =V + ^aV-^ {A + A^) + 0{a^) . (15) 

In this setting, it is still most natural to follow the notation of our Appendix 
and to define the operator of position of our (quasi)particle by the strictly 
diagonal matrix ([9]) acting in the "inaccessible" physical Hilbert space T-L^^^ 
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where the metric remains trivial, B'-^^ = I. Naturally, such a position oper- 
ator q may be equally well represented by the nondiagonal operator 



Q = Q-'qn 

acting either in the "accessible", unitarily equivalent physical Hilbert space 
"H^-^) (where, in our notation, Q = is Hermitian) or in the unphysical 
and unitarily non-equivalent auxiliary Hilbert space "H^^-* (where the same 
operator remains non-Hermitian of course, Q ^ Q^). 

At this point our approximate formula ( [T5|) for Q enables us to transform 
q into its image 



Q 



V-^aV-' {A + A^) +0{a^) 



V + ^aV~' {A + A^)+0{a^) 



VqV + [VqV~^ {A + A^) - V'^ {A + A^) qV] + 0{a^) 



which is, up to second-order corrections, a diagonally dominated tridiagonal 
matrix. The construction remains also feasible at A; > 1, provided only that 
we make use of a partitioning and replace the tridiagonality of Q^^^ by the 
block-tridiagonality of O^^'*. 

We may conclude that in our ad hoc physical Hilbert space 'H^'^^ assigned 
to the cryptohermitian Hamiltonian H^'^\a) the original basis \qs) [such that 
i\Qs))s' ~ ^ss'] lost its connection with the observable position. The role 
of the position-eigenstates is taken over by the kets \Xqs) representing our 
(quasi)particle(s) localized at a site s (with spatial coordinate g^). They must 
be constructed as eigevectors of our cryptohermitian operator of position, 

Q\xJ=qs\xJ, s = l,2,...,N. (16) 

This is a numerical problem, possibly simplified at small a when the operator 
of position predominantly couples the nearest-neighbor basis states. 

In the language of measurements the probability of finding the parti- 
cle (or quasi-particle) in question at the s— th site of the lattice may still 
be determined by formula ffTTj) . after its appropriate amendment of course. 



The form of this formula itself may even stay unchanged, provided only 
that we replace the old, Hermiticity-assuming definition of the wave function 
i/j{t,s) := {qs\i^{t)) (meaning the "probability-density quantity") by its new, 
cryptohermiticity-assuming update 

s+k 

^(e.)(t,5) = {x,.\Qk\m) = E {x,M<dk),A\m))s' ■ (17) 

s'=s—k 

With the time-independent parameter a 7^ a{t) this definition acquires, in 
the light of Eqs. ([6]) and ([7]), the final compact form 

s+k N 

^{e.)(t,5)= E (g.l(0/c),,,$^e^^"*L(5',a,K). (18) 

s'=s—k n=l 

This conclusion establishes the direct link between our cryptohermitian Hamil- 
tonian (jl]) and the measurement od the position. Formula f fTSj) indicates that 
the different choices of the physical metric will lead to the nonequivalent pre- 
dictions of the results of the experiment. Vice versa, the experiments may 
be used, in principle, as a valid source of information fitting the parameters 
in phenomenological metrics [7]. 

3 A few mathematical prerequisites 

On the level of the textbook quantum theory the main obstacle of accepting 
the input Hamiltonian matrix H^'^'^ 7^ (//*^^))^ (with real spectrum) lies in 
the technical nontriviality of making it compatible with the standard postu- 
lates of Quantum Mechanics. In Appendix we explained that this acceptance 
is based on the "hidden" Hermiticity 

//(^)(a) = {H^''\a)f (19) 

(a.k.a. cryptohermiticity [Sj). This property is defined in the less common 
but still entirely standard Hilbert space n^^l It may almost always be reread 
as the Dieudonne's [11] constraint f|T2|) imposed in the auxiliary, friendlier 
Hilbert space 'H^^\ 
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For our model (jl]) the major technical difficulties related to the necessary 
construction of the right-hand-side matrices in Eq. f|T9|) (i.e., of the eligible 
metrics) will be tractable due to the important observation that the metrics 
may systematically he constructed as special, (2/i;-|-l)— diagonal matrices (fT^ 
at any A; e {0, 1, . . . , - 1}. 

3.1 The exceptional local model at k = 

We shall assume that the parameter a and the dimension are both variable 
and, in principle, arbitrary while the number k of the nonvanishing diagonals 
in Eq. f lT3|) will remain, for technical as well as interpretation reasons, fixed. 
Naturally, the constructive analysis of the model should start from the trivial 
A; = 0. In this case it is well known that within the formalism described in 
Appendix, any tridiagonal complex matrix may be made compatible with 
the Dieudonne's Hermitization equation ^''^0 = HQ via the diagonal ansatz 



e 



11 



722 



7Af-l,Af-l 



6 



NN 



(20) 



for the metric with the positive matrix elements, 6jj > 0, j 
This observation leads to the following easy consequence. 



1,2,..., AT. 



Lemma 1. Hamiltonians H^^\a) of Eq. may be assigned the diagonal 
metrics Oo(a) of Eq. i^2U\) with normalization On = 1 and with the elementary 
matrix elements 



On 



atn = 2,3,...,N. 



(n-l)\ 



(a + 72 - 1) (a + n - 2) ... (a + 2) (a + 1) 



(21) 



Proof. In the Dieudonne's equation = /f© it is easy to evaluate the 
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left-hand-side product of matrices for any tridiagonal input 







02U22 


U 




(J 


U 




qfii 


(^2^22 




U 




n 
U 







C2O22 


03^33 


64^44 

















C*j^f_2^N~2,N-2 


fljV-l^Af-l,Af-l 






















(22) 



and to compare it with the right- hand- side tridiagonal matrix 



QHiN) 



Outti OiiCi 

^22^2 0220-2 ^22 C2 
6133 63 6*3303 






6*33 C3 






















^N-lM-lbN-l 0N-l,N-l(^N-l 0N~l,N-l^N~l 







e 



NN'JN 



Onno-n 
(23) 

The net result of this comparison is that the diagonal elements of the Hamil- 
tonian must be real and that the sequence of relations 



n+l,n+lbn+l 



r* 



n 



1,2, 



N-1 



(24) 



must be satisfied. For our present model 
leads to the proof of the Lemma. 



this observation immediately 

□ 



Remark 1. The simplicity of the above result made it useful in tests of a 
straightforward symbolic-manipulation algorithm which we intended to use 
for the computer-assisted solution of the Dieudonne's equation H'^Q = HQ 
using the various less elementary forms of the ansatz for {2k + 1)— diagonal 
metric Q = Qk- Identifying the present model with the /c = special case 
the algorithm produced the sequence of formulae 

6*11 n 2^11 



6'. 



22 



a + l 



e 



33 



(a + 2) (a + l) 
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and confirmed the validity of tlie above closed-form rigorous result fl2T]) and, 
ipso facto, also the reliability of the algorithm itself. 



3.2 The left eigenvectors of H^^\a) at k > 1 

As long as we have to work with the phenomenological model of quantum 
dynamics where H^^\a) ^ (^H^^\a)^ at all and a, it is insufficient to 
know just the a priori specified solutions (jTj) of the current time-independent 
eigenvalue problem ([5]). For multiple purposes it is also necessary to know the 
left eigenvectors of our Hamiltonian. Most often, we rather construct 
their duals defined as the right, usual eigenvectors of the conjugate, 
unusual Hamiltonian ^H^^\a)j , 

(i7W(a))^ 1^)) = )(a) 1^)) , n = 1, 2, . . . , A - 1 . (25) 

As long as we already know the energies, this is a simpler task. A rather 
subtle problem only emerges when we recollect that a suitable set of the left 
eigenbras of matrix H^^^a) has already been defined in Appendix 

(cf. Eq. (138!) ). This implies the necessity of the following proportionality 
relation 

l^r)):=0|^f^) = ie^)ena) (26) 

between the two alternative sets of the eigenvectors where the left-hand-side- 
state normalization is fixed while the right-hand-side-state normalization still 
remains open in Eq. (123]) . 

One of the consequences of the arbitrariness of the n— dependent complex 
constants Ki^\a) is particularly serious since the solutions of Eq. ([25]) 

are often used in the spectral-expansion definition of the metric [121 [TB| . 

Af-l N-1 

= E l^i''^)) ((^i^^^l = E l^'rT') Wn'^Ha)]' (e^l . (27) 

n=0 n=0 

This implies that up to an inessential overall multiplication factor (say, 
> 0), the metric (1271) contains, in general, an (A^ — 1)— plet of 
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free parameters K^J^\a) > 0, n = 1,2, . . . , N — 1. As a consequence, 
every given Hamiltonian H^^\a) may be assigned an (A^ — 1)— parametric 
multiplet of alternative metrics leading, generically, to independent physical 
predictions of the model. 

For this reason it is important to have, always, a sensible physics-based 
recipe for an efficient suppression of the latter ambiguity. In field theory one 
often requires the observability of the so called "charge" C for this purpose 
[5]. For quantum lattices using the nearest-neighbor interaction dynamics we 
proposed, very recently [10], another principle of suppression of the ambiguity 
of the metrics which is also being used here. This recipe is based on the 
requirement of the minimal (or at least "tunable" , controllable) non-locality 
of the model in question. In the latter setting, unfortunately, the direct use of 
spectral formula (127|) does not work at all. At the same time, the use of the 
metric-multiplication definition ( 126|) of the left eigenvectors li^n^^)) may be 
recommended as facilitated by the {2k + 1)— diagonal band-matrix structure 
of the metric, especially at the not too large k. 

4 The method of construction of the multi- 
diagonal metrics 

4.1 The sparse- matrix-expansion ansatz. 

Our preliminary computer-assisted experiments with solving Dieudonne's 
Eq. (fT2|) with Hamiltonian (jl]) revealed that its {2k + 1)— diagonal matrix 
solutions (fT3|) might be sought in a specific linear-superposition form 

e = ef V) = e^Z,^,,...,^Ja) = («) + a,vf \a) (28) 

with k G {1,2,...,A^ — 1} and with the diagonal metric of paragraph 13.11 
accompanied by a A;— plet of {2j + 1)— diagonal real and symmetric matrices 
Vj^\a). In fact, just the diagonal matrix Q'^\a) is a true, positive-definite 

14 



metric. Every other matrix "Pj (a) is only expected individually compatible 
with the Dieudonne's constraint, 

(^(A^)(^))t pW(^) = pW(^) ^(A^)(^) j = l,2,...,k. (29) 

In this overall framework the core of the exhaustive construction of the ex- 
pansion (128!) at a given non-locality k lies, obviously, in the constructions of 
all of the pseudometrics (!29|) with j < k and in having the real expansion 
parameters aj sufficiently small to keep the necessary positive definiteness of 
the diagonally-dominated sum (1251) unbroken. 

4.2 The extrapolation method of solving Eq. (1291). 

In what follows the systematic and explicit construction of pseudometrics 
Vj^\a) will be performed for all j < 3. This construction will proceed in 
three steps. 

In a preparatory step we select an integer value of subscript j > 1, insert 
a general (2j + 1)— diagonal-matrix ansatz for vf^\a) in Eq. fl29l) and, using 
a computer-assisted trial-and-error strategy, fine-tune an ansatz for V^^\a) 
in such a way that a maximum of its matrix elements not lying on its outer 
diagonals is being set equal to zero. 

In the second step of the algorithm we select a few N < Nq and, via the 
computer-assisted closed-form solution of Eq. (12^ we construct all of the 
corresponding matrix elements 9mm' = ^mSi^) concrete 
model this computer-assisted "experiment" revealed, after certain symbolic- 
manipulation factorizations and simplifications of the originally quite te- 
dious results, the prevailing A^— independence and sufficiently transparent 
and nicely factorizable a— dependence of the individual matrix elements of 
our separate pseudometric matrices Vj^\a), at the first few smallest j at 
least. 

In the third step the resulting sample of elements must be perceived and 
reanalyzed as a set which is generated by the comparatively elementary linear 
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algebraic set of equations fl29|) . These equations always possess, at any fixed 
A'^ and j, an implicitly recurrent structure. Such a reinterpretation of the 
algebra opens the possibility of skipping the intermediate step of the explicit 
determination of the individual recurrences as sampled, at j = 0, by Eq. fl2^ 
above. Especially at the higher values of j such an intermediate step did 
prove prohibitively complicated and, at the same time, redundant, especially 
due to the elementary nature of our Hamiltonian in question. 

Thus, in the third step we replace the tedious multidimensional recurrent 
generation of the closed algebraic formulae for the unknown multiplets of 
elements 6mm' by the extrapolation technique. In our model the latter recipe 
proved more efficient than the direct use of recurrences even at the lowest 
values of j. Naturally, having the closed form of extrapolation at our disposal 
at last, the final rigorous proof of its general compatibility with Eq. (!29|) by 
direct insertion remains straightforward. 

5 The model with minimal smearing (k = 1). 

In the general family of a— dependent tridiagonal metrics 



"numbered" by a real and not too large variable a we may assume, without 
any loss of generality, that the matrix v[^\a) is real and symmetric and 
that without loss of generality its element 6*11 may be chosen as vanishing. 



{a)+aVr\a) 



(30) 



9 



'12 



O12 022 ^23 



e. 



23 



6'. 



33 



'34 




6n-3,N-2 On-2,N-2 dN-2,N-l 

6n-2,N-1 dN-l,N-l 6n-1,N 

6n-i,n 6nn 



(31) 
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One of the most unexpected observations made during the brute-force con- 
struction of this pseudometric appeared to be the cutoff-independence of the 
result and, in particular, of the "last" element Ojsin 7^ 0. Thus the "initial" 
list of the elements 

^ - 2 . 8 

"22 — — T ) t'SS 



a + 1' (a + 2)(a + l)' 

^^^-~(a + 3) (a + 2)(a + l)' " ^ ' 

^ - 2 n 6 
f 23 — — —r ) t/34 



a + 1 ' (a + 2) (a + 1) 

obtained at < A^o = 4 proved sufficient for the extrapolation and for the 
final formulation of the general result. 

Lemma 2. Hamiltonians H^^\a) of Eq. may be assigned the tridiagonal 
family of metrics l[3U\) . The elementary matrix elements of pseudometric i^37\) 
are given by the cutoff-insensitive formulae 

Q ^ (^-1)! 

(a + n-l)(a + n-2)...(a + 2)(a+l) 

at n = 2, 3, . . . , A^, and 

n\ 

"nri+l 



(a n - 1) (a + n - 2) . . . (a + 2) (a + 1) 
at n = 1,2,..., N -I. 

Proof. In a parallel to the constructive proof of Lemma [1] the tridiagonal- 
metric-generated recurrent relations remain sufficiently transparent to admit 
the explicit identification and the straightforward proof of validity of their 
closed solutions by mathematical induction not only for our model (jl]) but for 
any tridiagonal input Hamiltonian ([H]). The details are left to the readers. □ 

5.1 The verification of extrapolation hypotheses 

The computation technique which we use here may be characterized as an 
interactive algorithm based on a systematic computer-assisted verifications 
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of the series of amended extrapolation hypotheses. It will remain applicable, 
mutatis mutandis, at any A; > 1, throughout our forthcoming constructions. 

The key technical problem will always lie in the determination of an 
appropriate ansatz for pseudometrics. In practice, the only way of finding 
the optimal ansatz seems to lie in a patient, brute- force solution of Eq. (fT2|) 
at a sequence of the smallest dimensions A^. That's what we will always have 
to do. Without the help of MAPLE, such a task would be rather difficult. 

In an illustrative verification of the result presented in Lemma [2] we may 
compare its predictions, say, with their two computer-generated counterparts 

24 192 
^^45 — 7 — , ON / — rrrTT — r~rT i *^55 



[a + 3) (a + 2) (a + 1) ' (a + 4) (a + 3) (a + 2) (a + 1) ' 

This comparison reconfirms the validity of Lemma [2l Moreover, these and 
similar "redundant" formulae may offer an insight and background for the 
estimates, say, of the numerical magnitude of the eigenvalues of the metric 
at larger parameters a and/or a. In the context of a different quantum- 
lattice model a confirmation of feasibility of such a nonperturbative search 
for the strong-coupling boundaries of the domain of positivity of the metric 
was mediated, e.g., by Table I of Ref. [H]. 

An additional merit of the interactive amendment of the tentatively ex- 
trapolated formulae via the a posteriori computer-assisted verification has 
been found in its speed. The initial tedious algorithm of the lengthy direct 
construction gets easily amended in the light of the extrapolation so that the 
calculations at the higher dimensions prove, paradoxically, quicker. Giving, 
in our illustrative example, the sequence of further elements 

120 



9. 



56 



(a + 4) (a + 3) (a + 2) (a + 1) 
1200 



66 



(a + 5) (a + 4) (a + 3) (a + 2) (a + 1) 
(etc) which further confirm the reliability of the formulae and open the way 
to the perceivably simplified proofs using direct insertions. 
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6 Pentadiagonal metrics {k = 2) 

In the general pentadiagonal metric of our model (jlj), 



(32) 



the only unknown matrix v'^\a) must again be real and symmetric, normal- 
ized, say, by the choice of ^(1,3) = 1 and finally, without loss of generality, 
admitting that the elements On and 612 vanish. 









^13 









^22 


6*23 






6*13 


^23 


^33 


^34 


6*35 



(^N-A,N-2 (^N-3,N-2 Qn-2,N-2 dN-2,N-l 0N-2,N 
0n-3,N-1 On-2,N~1 dN-l,N-l On-1,N 



On- 



N-2,N 



N 



t 



(N) 
NN 



(33) 



The hat-superscripted matrix element is exceptional and must be consid- 
ered manifestly cutoff-dependent. This has been revealed during the first-step 
calculations which gave the initial list of the cutoff-insensitive elements 

a + 2 ^ 4 (a + 5) ^ 18 (a + 8) 



722 



9: 



33 



^55 '- 
6*23 



a + 1 ' (a + 2) (a + 1 

96 (a + 11) 



6. 



(a + 4) (a + 3) (a + 2) (a + 1) 
4 „ 24 



a + 1 



6 



(a + 3) (a + 2) (a + 1) ' 

36 (a + 21) 
(a + 4) (a + 3) (a + 2) (a + 1) 
144 



34 



13 



1, 0; 



(a + 2) (a + 1) 
3 



ft 



45 



24 



(a + 3) (a + 2) (a + 1) 
12 



35 



a + 1' "^"^ (a + 2) (a + 1)' 
leading, by extrapolation and by its subsequent tests at a few > A^o = 5, 
to the following general result. 

Lemma 3. Hamiltonians H^^\a) of Eq. may be assigned the penta- 
diagonal family of metrics ^3^. The cutoff-insensitive matrix elements of 
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pseudometric are given by formula 

(n-l)(n-l)!(a + 3n-4) 



On 



(a + 72 - 1) (a + n - 2) ... (a + 2) (a + 1) 
at n = 2,3, . . . , N , and by 

2(n- l)n! 



Or, 



"""^^ (a + n - 1) (a + n - 2) . . . (a + 2) (a + 1) 
at n = 2,3, N — 1, and by 

Onn+2 



(n + l)\ 



2 (a + n - 1) (a + n - 2) . . . (a + 2) (a + 1) 
atn = 2,3,...,N -2. 

Proof. We leave the proof by direct insertion to the readers. 



□ 



Remark 2. The determination of the missing exceptional-element sequence 
t requires the use of a different approach, outlined in section [S] below. 



7 The metrics with seven diagonals, k = 3. 

In the seven- diagonal metrics 

ef )(a,a,/3,7) = G^^^la) + aP^la) + /^Pf ^(a) + 7^^^ (34) 

the missing real and symmetric component V^^\a) may be constructed along 
the same lines as above, starting from the assumptions On = O22 = 0, 
0{1, 2) = 0{1, 3) = and 0{1, 4) = 1 and from the heptadiagonal ansatz 



Pf'(a) 












Oii 












023 


024 









023 


033 


03A 


On-a,n-i 




Ou 


024 


03A 


Oaa 


0n-3,N~1 


On~3,N 




025 






On-2,N-1 


On-2,N 



■• 0n-3,N-1 0n~2,N-1 On_in-i 



N 







N-3,N 







N-2,N 



tjsi- 



f(N) 
''NN 



(35) 
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We now encounter the two specific, hat-superscripted matrix elements 
and i^^li jv which must be constructed in different manner (cf. section |H] 
below) . 

For the purposes of extrapolation we evaluated, this time, the following 
Nq = 6 set of the cutoff-insensitive matrix elements. 



a + 2)(a + l)' (a + 3) (a + 2) (a + 1) ' 

3a + 19 



^66 



^55 = -192 
= -4800 



(a + 4) (a + 3) (a + 2) (a + 1) 
a + 8 



(a + 5) (a + 4) (a + 3) (a + 2) (a+1) 



^45 = 36 
^56 = 240 



a + 1 ' (a + 2) (a + 1) 

a + 13 



(a + 3) (a + 2) (a + 1) ' 
a + 18 

a + 4) (a + 3) (a + 2) (a + 1) ' 

48 „ 360 



''^ a + 1 ' (a + 2) (a + 1) ' (a + 3) (a + 2) (a + 1) 

(714 — i , (725 — — , (736 



a + 1 ' (a + 2) (a + 1) 

leading, by extrapolation and its subsequent verification, to the following 

general result. 

Lemma 4. Hamiltonians H^^\a) of Eq. may he assigned the hepta- 
diagonal family of metrics [34\ )- The cutoff-insensitive matrix elements of 
pseudometric ^3^^ are given by formula 







2 (n - 1) {n 


-2) (n-1)! ( 


a + 5 n 


-6) 




3 (a + n- 1) 1 


[a + n -2)... 


(a + 2) 1 


(a + 1) 


at n = 2, 3, 


...,N, 


and by 












{n - 


- l)n\ [a + 5n 


-7) 






2 (a + n - 1) 


[a + n -2)... 


(a + 2) 


(a + 1) 


at n = 2, 3, 


...,iV 


— 1, and by 














(n - 1) (n + 1 


)! 






{a + n— 1)\ 


(a + n - 2) . . . 


(a + 2) 


(a+1) 
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at n = 2,3, N — 2, and by 

(n + 2)! 



Or, 



6 (a + n - 1) (a + n - 2) . . . (a + 2) (a + 1) 
atn = 2,3,...,N -3. 

Proof. Again, we leave the proof by direct insertion to the readers. □ 

8 The exceptional, cut-off-dependent matrix 
elements i^^^ 

8.1 The construction of tj\^j\^ cit k — 2. 

The mathematical nature of the determination of the exceptional elements in 
ansatz (133|) is special. Firstly, it is fairly clumsy to think about their values 
as resulting from recurrences since they vary with themselves. Secondly, 
the brute-force collection of data needed for extrapolation with respect to A^ 
is much more time-consuming. At the same time, one can extrapolate these 
data with respect to the integer variable A^ almost as easily as in the case of 
the other matrix elements. 

In such a situation we decided not to search for the rigorous proofs. Our 
readers will only be offered here the empirically multiply reconfirmed (i.e., 
for us, credible enough) results of the computer-assisted (though still rather 
time-consuming) extrapolations. 

Conjecture 1. At j = 2 and at all N = 3,4, . . . formula 
~(^) _ (A^-2)(A^-1)! (a + 5A^-4) 



2 (a + A^- l)(a + A^-2)...(a + 2)(a + l) 
defines the "last missing matrix element" of pseudometric (3^) . 

Remark 3. The validity of the latter formula (obtained by the extrapolation 

from the data at A^ < A^o = 5) has been made plausible, for us, by the direct 

MAPLE-mediated evaluation of the subsequent values up to 
~(9) _ 141120 (a -f 41) 

~ (a + 8) (a + 7) (a + 6) (a + 5) (a + 4) (a + 3) (a + 2) (a + 1) ' 
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8.2 The construction of tj\^-j\^ cit k — 3. 

The computer-assisted solution of Eq. fl2^ enabled us to evaluate the follow- 
ing set of results of symbolic manipulations, 

£(4) = _i6 " + '^ 

(a + 3) (a + 2) (a + 1) 



f(5) ^ 11a + 103 

(a + 4)(a + 3)(a + 2)(a + l) 



4? = -240 



= -960 



7a + 82 



(a + 5) (a + 4) (a + 3) (a + 2) (a + 1) 



17 a + 239 



''"^ (a + 6) (a + 5) (a + 4) (a + 3) (a + 2) (a + 1) 
which proved extrapolated as follows. 

Conjecture 2. At j = 3 and at all N = 3,4, . . . formula 

iiv) _ {N -3){N -1)\[{3N -A)a + 7 - IQN + 8] 



3{a + N -l){a + N -2) ...{a + 2){a + l) 
defines the diagonal missing matrix element of pseudometric 

Remark 4. Again, for us, the plausibility of this conjecture has been en- 
hanced by the MAPLE-mediated evaluation of the subsequent values up to 



= -80640 ■ 



23 a + 431 



(a + 8) (a + 7) (a + 6) (a + 5) (a + 4) (a + 3) (a + 2) (a + 1) 

A deeper meaning of the display of similar formulae may be seen not only 
in the misprint-control in the general formulae but also in the facilitated 
future identifications of the underlying recurrences if any. In addition, the 
transparency of the formulae (and, in particular, of their a— dependence) also 
underlines the formal simplicity of the metric itself, especially when we com- 
pare its matrix elements, say, with their complicated polynomial analogues 
as obtained in Ref. [TU] . 
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8.3 The construction 

In a way paralleling the preceding section we started from the data 

f(4) _ ^ a + 16 

''34 ~ ^ 



(a + 2) (a + 1) 
a + 23 



(a + 3) (a + 2) (a + 1) 



£(6) ^ ^^^^^ 

(a + 4)(a + 3)(a + 2)(a + l) 

''67 + 5) (a + 4) (a + 3) (a + 2) (a + 1) 

and formulated our last extrapolation hypothesis. 

Conjecture 3. At j = 3 and at all N = 3,4, . . . formula 
fiN) _ (iV-3)(iV-l)!(a + 7iV-12) 



'Af-lAf 



3{a + N-2){a + N-3)...{a + 2){a + l) 
defines the last missing off-diagonal matrix element of pseudometric ^3^). 

Remark 5. We may again display the two subsequent MAPLE-generated 
quantities, viz., 

4.\°) _ g4oo ■ 

(a + 6) (a + 5) (a + 4) (a + 3) (a + 2) (a + 1) 



and 



tla^ = 80640 



a + 51 



" (a + 7) (a + 6) (a + 5) (a + 4) (a + 3) (a + 2) (a + 1) 

The inspection of these formulae demonstrates that our model with the 
Laguerre-polynomial wave functions remains extrapolation-friendly even at 
k = 3. This is in contrast with the observations made in Ref. [H] where the 
complexity of the elements of metrics grew only too quickly between k = 2 
and k = 3. In loc. cit., we even failed to find a reasonable k = 3 formulae, 
due to a loss of any obvious guidance for the extrapolations. In the present 
model such an overall pattern seems to stay unchanged. Hence, one could 
also expect the reasonable feasibility of further, k > 4 constructions. 
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9 Summary and discussion 



A family of exactly solvable A^— site quantum lattices with a non-Hermitian 
nearest-neighbor interaction was proposed and studied. The energies ap- 
peared real so that each Hamiltonian has been made self-adjoint, i.e., stan- 
dard and physical in an ad hoc Hilbert space 7i^^^ where the inner product 
was defined via a metric G. The complete set of the eligible metrics has 
been shown numbered by a multiindex {k,ai, . . . , ak) in which the "degree 
of non-locality" k G {0, 1, . . . , A — 1} indicates that every 6 = Q(k,d) is a 
{2k + 1)— diagonal matrix. The other free parameters forming the A;— plets 
d must only be real and sufficiently small (otherwise, the metric could cease 
to be positive definite). In closed form the metrics were constructed and 
displayed for k = 0,k = l,k = 2 and k = 3. 

One of the main merits of the present quantum lattice model lies in a 
maximal suppression of its numerical aspects. First of all, the model does not 
need numerical methods for the solution of the time- independent Schrodinger 
equation. The reason is that the ket-eigenvectors of our special, solvable 
Hamiltonian H^^\a) were simply selected in advance (cf. Eqs. ([T]) or ([7])). 
We believe that from the point of view of fiexibility of the model such an 
apparently very strong a priori constraint has been more than compensated 
by the multiplicity of the eligible metrics &\!^^{ci) 7^ I- 

In a phenomenological context we explained that our model offers a new 
pattern of a "smearing" of the position of the lattice sites in the manner 
explained in paragraph 12.21 above. Such a form of nonlocality has been shown 
mediated by the metric. In accord with the "weakly nonlocal" interpretation 
of quantum systems as advocated in our recent paper [lOj we showed how 
the "site" of the lattice gets smeared over 2k + 1 neighbors, provided only 
that the range of smearing k is not too large. 

In practice people rarely select both the suitable Hamiltonian H = 
and the metric G = @(™p"*) as the two independent sources of informa- 
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tion about quantum dynamics. An important amendment of this limitation 
emerged in the hterature cca twelve years ago when a real boom of interest in 
quantum models with nontrivial metrics 7^ / has been initiated by Bender 
and his colleagues [T51 HSl [H] . These authors emphasized that the use of 
7^ / seems strongly motivated in quantum field theory. 

Incidentally, in a way explained by Mostafazadeh [9] and emphasized by 
Jones [IH] the latter studies were solely using, in our present language, the 
long-range metrics 0^^^ with the maximal possible subscript k = N — 1. This 
made the A; = — 1 recipe inapplicable in the unitary theory of scattering 
[20] . The core of the problem has been identified with the locality of 
the interaction f2T\. The restoration of the manifest unitarity of quantum 
scattering has only been achieved via the use of nonlocal interactions [22] . 

Fortunately, the difficulties of this type were successfully dealt with in 
various chain- interaction phenomenological models (cf. [23j). A number 
of constructive Q ^ I results appeared in the context of Bose-Hubbard 
[23] or Friedrichs-Fano-Anderso [25]) models. The appeal of these models 
in condensed-matter physics resulted in detailed descriptions of the tightly 
bound lattices of electrons [2S]), of the XXZ spin chains [27]) and, last but 
not least, of certain sophisticated experiments in optics [281 129] . 

In this context the distinctive feature of our present model lies in its 
exact solvability. Another relevant property of our present model is that its 
spectrum is safely real at all a > (cf. Eq. ([6])). We believe that such a 
spectrum might find useful applications, say, in a purely phenomenologically 
motivated fitting of measured energies. 

Although the similar idea has been also proposed and illustrated, via 
another model, in Ref. [H], the present spectrum (sampled, for reference 
purposes, in Table [T] above) looks much better suited for the fitting purposes. 
For various values of parameter a the present energy levels are regularly 
spread over a subinterval of the real axis which grows with the growth of 
the matrix dimension N. In contrast, the distribution of the bound-state 
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Table 2: A sample of the dimension- and parameter-dependence of the energy 
spectra {En{a)} for the quantum-lattice Hamiltonian of Ref. [H]. 



parameters 


energies 


A T 


a 


Eo[a) 


Ei[a) 




EN-2[a) 


EN-i{a) 


n 




i.U 


-.o44io554io 


-.3709690601 




.3709690601 


.6441855418 




2.0 


-.5622585222 


-.3463446402 




.3463446402 


.5622585222 




3.0 


-.5083600312 


-.3244642920 




.3244642920 


.5083600312 


9 


1.0 


-.6443127436 


-.3985980302 




.3985980302 


.6443127436 




2.0 


-.5626496595 


-.3794465124 




.3794465124 


.5626496595 




3.0 


-.5091239690 


-.3616352067 




.3616352067 


.5091239690 



energies of the model of Ref. |14j (sampled here by Table [2] for comparison) 
seems handicapped not only by its restriction to a fixed interval (—1,1) but, 
more seriously, by its not-well-motivated symmetry with respect to its center 
(shifted, conveniently, to the origin in Ref. [H]) and, even more seriously, by 
the weaker sensitivity of its extreme values to the changes of the dimension 
parameter A^. 

Besides the above-emphasized potential physical relevance of the exact 
solvability of our present family of Laguerrian bound-state models we would 
like to mention, in the conclusion, also the purely mathematical appeal and 
consequences of our closed-form constructions. 

Firstly, the unusual though still fully non-numerical solvability of our 
model could serve, in principle at least, in the role of the initial zeroth ap- 
proximation in perturbation theory. After all, just a very few exactly solvable 
models with metrics G 7^ / exist on the market. More attention has only been 
paid to a few exceptionally simple quantum lattices with certain extremely 
elementary point-like interactions [30| [31]. 

Secondly, one should emphasize that it was certainly unexpected that 
the innocent-looking freedom in the choice of the proportionality constants 
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Kn (a) in relation ( 126|) proved able to render the general definition of met- 
ric ( 1271) compatible with the fairly strong simplification requirement of the 
sparse- matrix structure of its {2k + 1)— diagonal representations Q\!^^{a) as 
specified by Eqs. (113]) and (EH]). 

From the point of view of the theory of classical orthogonal polynomials 
the latter result implies that every set of column vectors ([7]) of Laguerre 
polynomials may be assigned many different biorthogonalized sets of row 
vectors. They may be formed of the "ketket" eigenvectors \tp)) of [//*^^)(a)]^ 
and classified by the same multiindex {k,d} as the metrics, therefore. Indeed, 
we have 

lV^fH«))) = 0S(«)lV^fna)) 
or, in the componentwise notation of section 12. 2^ 

mm{N,s+k) 

(ieHa)))),= E (e2(a))^^,(ieH«))). (36) 

s'=max(l,s— fc) 

where s = 1, 2, . . . , and where matrices ©^^''(a) are now at our disposal in 
closed form fl25]) . for < A; < 3 at least. 

Thirdly, the variability of the /c— parametric {2k + 1)— diagonal metrics 
Q\f^]{a) could prove useful during the model-building in which the energies 
are prescribed by our "Laguerrean" input Hamiltonian H^^\a) but in which 

several other matrices of observables might be required selfadjoint in at least 

(s) 

one of the eligible Hilbert spaces of states T-Ll typically, via the fitting of 
the k free parameters a. 

Fourthly, on the level of methods the present approach to solvability 
could open new ways of circumventing the difficulties which emerged during 
the study of scattering with H ^ [HI [19]. The use of the non-numerical, 
exactly solvable models with tridiagonal Hamiltonians and large N oo 
seems to have been quite effective in this context [231 132] . 

Fifthly, the quick progress achieved during the study of PT— symmetric 
solvable differential equations [53] has not been followed by the sufficiently 
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rapid progress in understanding of the related metrics. One of the reasons 
hes in the virtually prohibitive technical obstacles |9]. In contrast, the recent 
transition to PT— symmetric difference Schrodinger equations has been ac- 
companied by the comparatively quick success in finding the comparatively 
extensive sets of metrics G compatible with a given H [lOl [321 

Sixthly, a deeper understanding of the problem of ambiguity of the as- 
signment of a metric to an input Hamiltonian has been reached here via the 
matrix Schrodinger equations. We pointed out that the old paradoxes [7] 
are finding new resolutions, say, due to the feasibility of a "smearing" of the 
coordinates at < oo. Due to the variability of k in our model we showed 
how the smeared localization could be studied by the experimental measure- 
ments, in principle at least (cf. also the possible application of this idea to 
cryptohermitian guantum graphs as mentioned in Refs. [Mf [35]). 

Last but not least, our continuing attention paid to the exactly solvable 
models could clarify the possibilities of the practical use of the concept of the 
hidden Hermiticity (cryptohermiticity) also in connection with the unitary 
time evolution where the time-dependent metrics emerge. As long as this idea 
leads to many new technical challenges [3S] , our present non- numerical model 
could re-demonstrate, in the nearest future, its relevance and importance also 
in this new context. 
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Appendix: Three Hilbert spaces 1-L'^-^\ and 
Ti^^^ and the concept of hidden Hermiticity 

Whenever a non-Hermitian Hamiltonian matrix H ^ H'^ with real spectrum 
is declared physical, the Hilbert space of states cannot be represented by the 
most common A^— dimensional version of The latter space may only play 
an auxiliary role (this space without any immediate physical interpretation 
will be denoted here by the symbol "H*-^-' where the superscript ^^'^ may stand 
for "friendly" as well as for "false" [5]). 

For the work with Hamiltonians H ^ there exists a trick introduced in 
physics, presumably, by Scholtz et al |7j. Its essence lies in the replacement 
of inappropriate "H^^-* by a unitarily non-equivalent, "standardized" Hilbert 
space H^^K By construction, the two Hilbert spaces share the same set of 
kets (in our notation, \ip) G V*-"^-* = V*-^-* = V). The only distinguishing 
feature is that the usual Hermitian conjugation, i.e., the most conventional 
transition to the duals, T^^^ : V — > V (mapping 

r(^): 1^) ^ (^1 in •H(^) (37) 

in the current Dirac's notation) must be replaced, in H^^^ by the generalized 
conjug ation r^"^) : V (= V^^^ ) ^ [Vf^\ The latter mapping generates the 
different duals defined in terms of a suitable matrix G of metric [8\, 

An instructive illustration of this scenario may be found in paper [7] where 
the kets |^) e V = V^^) = V^^) were just bosonic, imperfect representatives 
of the true nuclear (i.e., fermionic) state vectors 

=n\i;) G W (39) 

which lied in the textbook Hilbert space n^^^ with metric 0^^^ = I. The 
so called Dyson's map 1) : V — ?■ W of bosons upon fermions f p9|) was chosen 
non- unitary so that ^ipa\'4'b>- 7^ (^al'^b) since ^ Qr^ . 
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In opposite direction the non-unitarity of the map f2 enables us to require 
that 

-<A\ll^b^ = = ((^al^b) . (40) 

This formula connects the Dyson's map and the metric by the most important 
relation 

e = Q^VL. (41) 

The same picture of reality is obtained in both of the alternative physical 
Hilbert spaces of states 'H*-'^^ and l-i^^^ which are, by construction, unitarily 
equivalent. 

In the light of definition fl5^ the result of action of H upon |a), i.e., a new 
vector \h) = H\a) appears in equivalent relation VL^^\h>- = HQ^^\a>- , i.e., 
we have \b>- = [}\a>- where we abbreviated f) = QHQ~^. The latter image 
of Hamiltonian must be self-adjoint in 'H^^\ i.e., the operator [) = VlH^l^^^ 
must be equal to its conjugate in 'H''^\ viz., to [)''' = (il^^)^ H'^il^ . Such a 
constraint imposed in H^^^ is strictly equivalent to formula ([T2|) valid in the 
other two spaces H^^^ and H^^K 

For the Hamiltonian operator we have to distinguish between its apparent 
non-Hermiticity H W m. unphysical space "H*-^-* and the true, "hidden" 
Hermiticity in the physical Hilbert space "H*^"^-' a.k.a. cryptohermiticity. In 
the latter space the Hermitian conjugate of H (let us denote it by the symbol 
H^) is defined, consistently, by the prescription 

:= Q-^H^ e in U^^^ . (42) 

In this language, Eq. fll2p is precisely the disguised condition H = of 
the Hermiticity of the Hamiltonian at a fixed G. This scenario has been 
used, e.g., in Ref. [7j. On the contrary, whenever we select the concrete form 
of matrix H in advance, relations ([T^ must be reread as the Dieudonne's 
(incomplete) set of linear equations for the matrix elements of = <d{H). 
This is precisely the approach used in our present paper. 
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